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1 Abstract 

New preliminary BaBar results for rare leptonic decays K~vV and —>■ 

are reported. Using data collected at the T(4S') with the BaBar detector, no evidence 
for a signal was found yielding the corresponding upper limits at the 90% confidence 
level: B{B~ K-yV) < 9.4 x 10"^ for 50.7 fb"^ , B{B^ e+e") < 3.3 x 10"^ 
^ < 2.0 X 10-^ and B{B^ e^/j,^) < 2.1 x 10"^ using 54.4 fb~^ . 

2 B- ^ K-pv 

The exclusive decay B~ K~vV is characterized by the absence of the long distance 
contributions and by the fact that the effective Hamiltonian is represented in the 
Standard Model (SM) by only one operator. In the SM the decay proceeds via the 
W box and Z penguin diagrams as can be seen in fig |1|. This mode probes the quark 
mixing parameter The SM prediction is B{B~ R-uV) = S.StH x 10"^ 
and the previous best hmit obtained by CLEO is B{B~ K^vV) < 2.4 x 10"'' [|]. 
Enhancement beyond the SM can arise from various types of models . The highly 
constrained models where the existing bounds on other Flavor Changing Neutral 
Current processes imply that the rate for B~ K~vV can not exceed the SM 
prediction by more than a factor of a few. In this category we have for example the 
Multi Higgs Doublet Models and the Left Right symmetric models. Enhancement by 
one or two orders of magnitude can arise from models with an extra vector-like down 
quark or models with leptophobic Z bosons 0. Last are the unconstrained models 
where the couplings responsible of enhancing B~ K~vV are to a large extent 
independent of current existing experiment bounds. An example which belongs to 
this category is supersymmetric models without R-parity. 

The presence of two neutrinos in the final state precludes the use of any kinematic 
constraints on the signal B meson. The strategy adopted for the BaBar analysis 
was to look for a charged kaon with momentum p* > 1.5 GeV in the T(4S) frame 



1 




B" K" B" K" B" K" 

Figure 1: Standard Model Feynman diagrams for B~ K~vV. 



recoiling against a semileptonic decay D l^u{X) where X represent either a 

null or a photon or a vr*^ from higher mass charm states. The exclusively reconstructed 
B meson is referred in the following as the tag B. The BaBar standard practice of 
blind analysis was followed to prevent biases. The low multiplicity of the signal 
decay reduces the combinatorieal background in the tag B reconstruction allowing 
the semileptonic decay B^ 'lf'i^v{X) to be cleanly reconstructed. The is 
reconstructed in the K^ti~ , K'^h^ti^tt'^ and K'^ir'ir'^ modes. This method results in 
roughly 0.5% of tag reconstruction efficiency. The data used in the analysis consist 
of 50.7 fb"^ collected at the T(4S) resonance corresponding to 56.3 x 10^ BB events 
and 6.4 fb~ collected just below BB threshold. The kinematics for the B~ K^uV 
decays in the simulation is based on the form factor model in 

Hadronic events were selected once an electron or muon with a momentum above 
1.3 GeV in the T(4S') rest frame was identified. Than D candidates were recon- 
structed in one of the K^tt' , K^ti'ti'tt^ and K^tt'tt^ decay modes. The kinematic 
requirement on the angle between the B and the reconstructed , calculated in 
the T(4S) frame was used to suppress background and restrict the kinematics of the 
!D^£"'" to be consistent with coming from a semileptonic B decay. The requirment 
—2.5 < cos Di < 1-1 was imposed, using 

cosObdi = — a n\ 



where Eb and \pb\ are respectively the energy and magnitude of the momentum of 
the B meson in the T(4S) frame. 

A special double tag sample was used to extract a correction to the efficiency 
calculated from signal Monte Carlo (MC). The sample was reconstructed by first 
finding a suitable D^i~^ candidate where the decays to K~^7i~, and then looking 
for a second Di candidate in any of the accepted D modes. The observed rate 
of double tags per fb~^ in the data is 0.85 ± 0.11 times the rate in the simulation 
leading to a correction to the the signal efficiency by a factor 0.92 ± 0.06 where the 
uncertainty is taken as a systematic error. The uncertainty in the efficiency of several 
of the selection criteria were also studied using the double-tagged sample. The total 
relative uncertainty on the selection efficiency was found to be 5e/e = 8.7 % where 
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the tagging efficiency and (-E'lcft) the remaing neutral energy after the tag B and its 
daughter were removed contribute the most. 

The distribution of events in the search plane defined by the variables[] Emi and 
(m^) — m^)/a^ is shown in fig We observe 2 events in the signal box, defined by 
the requirements i^ieft < 0.5 GeV and {m^ — m^) < Serf}. The expected background 
from the MC is 2.2 events. The background at present appears to be mostly cc events. 

The Poisson upper limit calculated at 90% C.L. without background subtraction 

is 

B{B~ ^ K-i^V) <9A-10-^ (2) 
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Figure 2: The distribution of events in the {{mu — mj^)/a^ -Eieft) plane for on-peak 
data, generic BB and continuum MC and signal MC. In the generic MC plot the blue 
circles show the contribution from BB events, the green squares show the contribution 
from cc and the red triangles show the contribution from uu/dd/ss. The generic MC 
contributions needs to be scaled by a factor of 1.09, 2.21 and 3.56 for the BB, cc and 
uu/dd/ss contributions respectively to correspond to the on-peak data luminosity. 



3 ^ l+l- 

B^ l^l^ proceeds in the SM through the three dominant W box and Z penguin 
diagrams shown in fig. ^ Even though these diagrams are similar to the one that lead 
to B~ — > K^uV these decays are further helicity suppressed by factors of mf. Both 
B^ l^l^ and B^ K^vV FCNC transitions provide an essential opportunity 
to test the SM and offer a complementary strategy in the search for new physics 
by probing the indirect effects of new particles and interactions. The SM theoretical 
branching ratio predictions are 1.9 x 10~^^ for B^ — > e"'"e~, 8.0 x 10~^^ for 5° — > 
and null for B^ — ^ e^fi"^ although recent neutrino mixing experiment results suggest 
that the branching ratio would be less than 10~^^. To date these decays have not 

^Thc quantities and af} are the mean and sigma from Gaussian fits to the invariant mass 
spectrum. Separate values are calculated for each decay mode in data and simulation. E'loft is the 
remaing neutral energy after the tag B and its daughter were removed 
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Figure 3: Standard Model Feynman diagrams for 




been observed and the current best limits from CLEO and Belle are summarized 
in table |^ 0-0 • Since these processes are highly suppressed in the SM they are 
potentially sensitive probes of physics beyond the SM. Models such as MHDM with 
Natural Flavor Conservation and large tanf3 can give up to an order of magnitude 
enhancement The enhancment in models with an extra vector-like down quark 
can be to up to two orders 0. Large enhancement may also arise from Minimal 
Supersymmetric Models with large tanf] or Supersymmetric models without R- 
parity. 

The decay offers a very clean experimental signature in the BaBar 

detector. The two high momentum leptons can be measured precisely and identified 
with high purity in the detector. Relatively low backgrounds arise from the contin- 
uum consisting mostly of non-resonant e^e~ —>■ qq production where q = u,d, s, c. 
The main contribution in the case of B^ — > e~^e~ is from pairs of real electrons in 
cc production; the contribution of misidentified hadron-electron pairs is negligible. 
Two-photon events contribute a significant background. Misidentification of muons 
is significantly more important in the case of B evidenced by an in- 

creased background expectation from uds events. For the B^ ^'^fi~ channel, the 
background from two-photon processes is negligible. 

The main measurement criteria are used to suppress different background events. 
The multiplicity cut Nmuit = Nxrk + > 6,E.y > 80MeV suppresses radiative 

Bhabha events while maintaining a higher efficiency than a simple stringent multi- 
plicity cut. The tracks are restricted to the central region of the detector using a polar 
angle cut which is efficient in removing QED background due to its strong depen- 
dence on the polar angle. As mentioned above, the main background process comes 
from continuum events which exhibit a two jet structure and produce high momen- 
tum approximately back to back tracks satisfying the requirements imposed on our 
candidate events. Two different shape variables suppress this kind of background: 
I cos I < 0.84 where 6t is the angle between the thrust axes of the 5-candidate 
and the rest of the event and the event thrust magnitude |T| < 0.9. The lepton 
pair is selected by simultaneous requirements on the energy difference AE and the 
energy-substituted mass rriEs defined in the following. The invariant energy differ- 
ence of the i?-meson candidate and the energy-substituted mass ruEs are calculated as 
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AE 



PB-Pi- s/2 



niES 



./2 + p^-p,)VK2 



(3) 



where ^/s is the CMS energy, ps and pi denote the four- momenta of the 5-meson 
candidate and the initial state. 

The signal box is defined in the {mEs AE) plane and was obtained for each of 
the mode separately using an upper limit optimization. The size of the 

signal and Grand Sideband (GSB) was chosen to be roughly [+2, —2]a of the expected 
resolution in AE and [+2, —2\a for niEs- In the cases of e~^e~ and B^ e^yu"^, 

the signal box size in AE were relaxed to roughly [+2, —3]a and [+2, — 2.5]cr to 
account for increased amounts of final state radiation and bremsstrahlung. 

The resulting efficiencies for the three B^ I'^l" are given in table |l]. The 
number of events observed in the GSB appears in figure ^. In order to estimate the 
number of background events in the signal box an unbinned maximum likelihood fit 
was done using an Argus fit for the mEs distribution and an exponential fit for AE 
distribution. The results are given in table |I[ There are three sources of systematic 
uncertainties: The normalization, the signal efficiency and the background estimate. 
The systematic uncertainties for the different variables were estimated using a control 
sample oi B ^ J/'^K^ with J/^ ^ . Since there is no control sample for 
B^ — > e^yU^ the error for B^ e^e~ was assigned to B^ e'^jj^ and the total 
systematic errors amount to 8.6%, 5.2% and 8.6% for B^ — > e+e~,-B° — > /i"*"/!" 
B^ — i> e^/i^ respectively where the main contribution is from niEs and AE. 
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Figure 4: Unblinded {mEs, AE) distributions with 25, 26 and 37 events in the GSB 
for B^ e~^e~ ,B^ — > fi'^fi~ and B^ e^fi^ respectively. 

After unblinding, we observed 1 event in the signal box for B^ — > e~^e~ and none 
for the rest of the channels are summarized in table |l]. The unblinded {rriEs, AE) 
distributions for the three channels are shown in figs [Q. The observation are com- 
patible with the expected background. We do not perform background subtraction 
for the determination of the branching fraction (upper limit). 

The upper limits on the branching ratios for B^ l^l^ obtained at the 90% 
confidence level are summarized in table 
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Channel 


TV 

' exp 


Nobs 


Nbg 


£[%] 


UL (90% CL) 


e+e^ 


1 X 10-^ 


1 


0.60 ± 0.24 


19.3 ± OAOstat ± l.QOsyst 


3.3 X 10"'^ 


B^ 


4 X 10-3 





0.49 ±0.19 


18.8 ± 0.28,tat ± 2msyst 


2.0 X lO"'^ 









0.51 ±0.17 


18.3 ± 0.38,jai ± l.^Osyst 


2.1 X 10"^ 



Table 1: Summary of analysis. Nexp is the number of expected signal events assuming 
a branching fraction of 10"^^. Nobs is the number of observed events in the signal 
box. Nbg is the expected number of background events in the signal box. 



4 Results 

A summary of the BaBar preliminary results for the rare leptonic B decays are given 
in table ^ in comparison with CLEO and Belle results. 



Mode 


CLEO 


Belle 


Babar 


B{B- K-uV) 


2.4 X 10-^ 




9.4 X 10"^ 


BiB'^ e+e-) 
b\b^ e+i2+) 


8.3 X 10-^ 
6.1 X 10"^ 
15.0 X 10"^ 


6.3 X 10"^ 
2.8 X 10-^ 

9.4 X 10"^ 


3.3 X 10-^ 

2.0 X 10'^ 

2.1 X 10"^ 


Luminosity 


9.1 fb"' 


21.3 fb~' 


54.4 fb"' 



Table 2: Summary of BaBar, CLEO and Belle results for rare leptonic B decay. 
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